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The eigenvalues of the Dirac operator on a curved spacetime are diffeomorphism-invariant functions 
of the geometry. They form an infinite set of "observables" for general relativity. Recent work of 
Chamseddine and Connes suggests that they can be taken as variables for an invariant description 
of the gravitational field's dynamics. We compute the Poisson brackets of these eigenvalues and find 
them in terms of the energy-momentum of the eigenspinors and the propagator of the linearized 
Einstein equations. We show that the eigenspinors' energy-momentum is the Jacobian matrix of 
the change of coordinates from metric to eigenvalues. We also consider a minor modification of the 
spectral action, which eliminates the disturbing huge cosmological term and derive its equations 
of motion. These are satisfied if the energy momentum of the trans Planckian eigenspinors scale 
linearly with the eigenvalue; we argue that this requirement approximates the Einstein equations. 
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I. INTRODUCTION 

One of the important lessons that we learn from gen- 
eral relativity is that fundamental physics is invariant 
under diffeomorphisms: there is no fixed nondynamical 
structure with respect to which location or motion could 
be defined. A fully diffeomorphism-invariant description 
of the geometry has consequently long been searched in 
general relativity; but so far without much success. As 
emphasized by many authors, such a description would 
be precious for quantum gravity 0| . 

Recently, Alain Connes' intriguing attempt of using 
the particle physics standard model for unraveling a mi- 
croscopic noncommutative structure of spacetime 
has generated -in a sense as a side product- the remark- 
able idea that the curved-spacetime Dirac operator D 
codes the full information about spacetime geometry in 
a way that can be used for describing the dynamics of 
the latter. Indeed, not only the geometry can be recon- 
structed from the (normed) algebra generated by (the 
inverse of) D and the smooth functions on the manifold, 
but the Einstein-Hilbert action itself is approximated by 
the trace of a simple function of D ||^. In its simplest 
reading, this result suggests the possibility of taking the 
eigenvalues A" of the Dirac operator as "dynamical vari- 
ables" for general relativity. Indeed, these form an infi- 
nite family of fully four-dimensional diffeomorphism in- 
variant observables, precisely the kind of object that was 
long searched in relativity. This approach has limita- 
tions. The most serious of these are that so far it works 
for Euclidean general relativity only (see ||^ for some at- 
tempts to overcome this problem), and somewhat unre- 
alistic predictions for the bare couplings [Q. However, it 
definitely opens a new window on the study of the dy- 
namics of spacetime. 



In order to use these ideas in the classical or in the 
quantum theory, one must translate structures from the 
metric variables to the A" variables. In particular, one 
needs information on the constraints that the A" 's satisfy 
if they correspond to a smooth geometry, and on their 
Poisson brackets. The difficulty is that the dependence 
of the A"'s on the geometry is defined in a very implicit 
manner, and these quantities seem too hard to compute. 

In this letter, we show that these difficulties can be 
circumvented. Following some earlier results in (valid 
only for the non 4d-invariant eigenvalues of the fixed-time 
Weyl operator), we derive here an expression for the Pois- 
son brackets of the Dirac eigenvalues. Rather surpris- 
ingly, this expression turns out to be given in terms of 
the energy-momentum tensors of the Dirac eigenspinors. 
These tensors form the matrix elements of the Jacobian 
matrix of the change of coordinates between metric and 
eigenvalues. The Poisson brackets are quadratic in these 
tensors, with a kernel given by the propagator of the lin- 
earized Einstein equations. Thus, the energy-momentum 
tensors of the Dirac eigenspinors turn out to be the key 
tool for analyzing the representation of spacetime geom- 
etry in terms of Dirac eigenvalues. 

We study also the Chamseddine- Connes spectral ac- 
tion. In the form presented in |^] this is quite unrealistic 
as a pure gravity action, because of a huge cosmological 
constant term that implies that the geometries for which 
the action approximates the Einstein-Hilbert action are 
not solutions of the theory. However, a very small mod- 
ification of the action eliminates the cosmological con- 
stant term. We derive the equations of motion directly 
from the (modified) spectral action. We argue that they 
amount to the requirement that the energy momenta of 
the high mass eigenspinors scale linearly with the mass, 
and that this requirement approximates the vacuum Ein- 
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stein equations. 

These results suggest that -even independently from 
its application to the standard model- the Chamseddine- 
Connes gravitational theory can be viewed as a manage- 
able gravitational theory by itself (see also ) , possi- 
bly with powerful applications to classical and quantum 
gravity. It reproduces general relativity at low energies; it 
is formulated in terms of fully difFeomorphism invariant 
variables; and, of course, it prompts fascinating exten- 
sions of the very notion of geometry. 

II. GR IN TERMS OF EIGENVALUES 

Consider Euclidean general relativity (GR) on a com- 
pact 4d (spin-) manifold without boundary. We work 
in terms of the tetrad field Ej^{x). Here /i = 1 . . .4 are 
spacetime indices and I = 1 ... A are internal Euclidean 
indices, raised and lowered by the Euclidean metric Sjj. 
The metric field is given by gfj,u{x) — E'^{x)E^ j{x), 
and is used to raise and lower spacetime indices. The 
spin connection w^j is defined by the equation d[fj,E^^ — 

"^^^ dynamics is defined by the action S[E] — 
J d'^Xy^R, where g and R are the determinant and 
the Ricci scalar of the metric. 

In spite of a curiously widespread popular belief of 
the contrary, phase space is a covariant notion: the co- 
variant definition of the phase space is as the space of 
the solutions of the equations of motion, modulo gauge 
transformations ||ll|. In the theory considered, the gauge 
transformations are given by 4d diffeomorphisms and by 
the local internal rotations of the tetrad field. Thus, the 
phase space T of GR is the space of the tetrad fields E 
that solve the equation of motion (Einstein equations), 
modulo internal rotations and diffeomorphisms. F can 
be identified with the space of the gauge orbits on the 
constraint surface and with the space of the Ricci flat 
4-geometries. 

We shall use the following notation. We denote the 
space of smooth tetrad fields as £; the space of the orbits 
of the gauge groups -diffeomorphisms and local orthogo- 
nal tetrad rotations- in £ as (these are "4-geometries" ); 
and the space of the orbits that satisfy the Einstein equa- 
tion as r (these are the Ricci-fiat 4-gcometries, which 
form the phase space of GR) . We call functions on F "ob- 
servables". Observables correspond to functions on the 
constraint surface that commute with all the constraints. 

We now define an infinite family of such observables. 
Consider spinor fields ip on the manifold and the curved 
Dirac operator 

D^ij'E^;{d^ + u;^jKl'l''), (1) 

which acts on them. Here 7 are the (Euclidean) hermi- 
tian Dirac matrices. For each given field E, the Dirac 
operator is a self-adjoint operator on the Hilbert space of 
spinor fields with scalar product 



(V',<^)- j d'x^ il:{x)c^{x). (2) 

where the bar indicates complex conjugation, and the 
scalar product in spinor space is the natural one in C*. 
Therefore, D admits a complete set of real eigenvalues 
and eigenfunctions ("eigenspinors"). Since the manifold 
is compact, the spectrum is discrete. We write 

I^V;" = A"?A". (3) 

Here and below, n = 0,1,2... is an index, not an expo- 
nent. We convene to label the eigenvalues so that A" 
is non decreasing in n, namely A" < A""*"^ (each eigen- 
value is repeated according to its multiplicity). In order 
to emphasize the dependence of Dirac operator, eigenval- 
ues and eigenspinors on the tetrad field, we use also the 
notation 

Z?[S] = A"[^] (4) 

where the dependence on the tetrad is indicated explic- 
itly. 

A" [E] defines a discrete family of real functions on the 
space £ of the tetrad fields 

A" : £ — ^R 

E^X'[E]. (5) 

Equivalently, they define a function A from £ into the 
space of infinite sequences R°° 

A : £ — >R°° 

E^{X\E]}. (6) 

Since we have chosen to order the A" 's in non-decreasing 
order, the image of £ under this map, which we denote 
as X{£) is entirely contained in the cone A" < A""''^ of 

Now, for every n, the function A" is invariant under 
4d diffeomorphisms and under internal rotations of the 
tetrad field. Therefore the functions A" are well defined 
functions on Q. In particular, they are well defined on F: 
they are observables of GR. 

Two metric fields with the same A"'s are called 
"isospectral" . Isometric (that is, gauge equivalent) E 
fields are isospectral, but the converse might fail to be 
true flj,^. Therefore A may not be injective even if re- 
stricted to Q. The A"'s may fail to coordinatize Q. They 
may also fail to coordinatize F. However, they presum- 
ably "almost do it" . In the following, we explore the idea 
of analyzing GR in terms of the set of observables A" [E] . 

III. SYMPLECTIC STRUCTURE 

The phase space F is a symplectic manifold and a Pois- 
son brackets structure is defined on the physical observ- 
ables. We now study the Poisson brackets {A", A™}. 
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To this purpose, we first construct the symplectic 
structure on T. This can be written in covariant form 
following Ref H. A vector field X on r can be written 
as a differential operator 



definition of the propagator of the linearized Einstein 
equations over the background E (in the chosen gauge). 
For instance, let us chose the surface E as x"* = and fix 
the gauge with 



J d'xX'^ix)[E] 



(7) 



X' 



1, X',=Q, Xl^l, Xl^Q. 



(14) 



where Xl^{x)[E] is any solution of the Einstein equations 
linearized over the background E. The symplectic two- 
form n of GR is given by 



where a = 1,2,3 and i = 1,2,3. Then equation ( |T^ ) 
becomes 

F^{x, t)= [ d'y (P:5(x, t; y, 0)foFi{g 0)), (15) 



n{x, Y) 



— - — / d'^a np[X^ Y^Wuv^^'''^" (8) where we have used the notation x — {x^.x'^.x^) and 
327rG t — x'^. and the Drooaeator can be easilv recomized. 



where [X^ Yj] = [X^ V, Yj ~ F/ V, Xl\, from 
now on we put 327rG = 1. Both sides of ^ are functions 
of E, namely scalar functions on F; this E is used to 
transform internal indices into spacetime indices. Here 
S : a I — > x{(t) is an arbitrary three-dimensional "ADM" 
surface, and Up is the normal one-form to this surface. 
The coefficients of the symplectic form can be written as 

nZ{x,y)= f np[d{x,x{a))^rS{y,x{a))]e]j^e-'"^\ 

(9) 

The Poisson bracket between two functions / and g on 
F is given by 



{f.g} 



d^x I d'y Pllix^y) 



Sf 



Sg 



6EUX) SEUyY 



(10) 



where P^{^(a;, y) is the inverse of the symplectic form ma- 
trix. It is defined by 



d^y Pli{x, y) ny^iy, z) = <5(x, z) Sj,. (11) 



Since the symplectic form is degenerate on the space of 
the fields (it is non-degenerate only when restricted to the 
space of equivalent classes of gauge-equivalent fields) , we 
can only invert it on this space by fixing a gauge. Let us 
assume this has been done. More precisely, integrating 
the last equation against a vector field F^{z) that satis- 
fies the linearized Einstein equations over i?, we have 



Jd^yjd^z P^J{x,y) ny^iy^z) F^ [z) ^ 
^Jd''zS{x,z)SPSi,Fp^iz), 



(12) 



Integrating over the delta functions, and using (|^), we 
have 

J^d^a np[Pli{x,x{a))^ pF^{x{a))yj^^e^'^^^ = 

- Fl{x). (13) 

This equation, where F is any solution of the linearized 
equations, defines P. But this equation is precisely the 



t — x'^, and the propagator can be easily recognized. 

Next, we need the functional derivative of the eigen- 
values with respect to the metric. The variation of A" for 
a variation of E can be computed using standard tech- 
niques for first order variations of eigenvalues; (standard 
time-independent quantum-mechanics perturbation the- 
ory). For a self-adjoint operator D depending on a pa- 
rameter V and whose eigenvalues are nondegenerate, we 
have 



i7A" d 



(16) 



In our situation, for generic metrics with nondegenerate 
eigenvalues we have that 



A A" 



'SEf^ixY 



(17) 



SEUx) 



SElix) 
5 

seJJ^) 



5El{2 



s 



5El{x) 
T"^(x). 



(18) 



Now, T"j{x) is nothing but the usual energy- momentum 
tensor of the spinor field ■0" in tetrad notation (see for in- 
stance [0). Indeed, the usual Dirac energy- momentum 
tensor is defined by 



T^ix) 



SE^,{x 



(19) 



where ^Dirac = / ^/g {ipDip — X'ipip) is the usual Dirac 
action of a spinor with mass A. 

We have shown that the energy-momentum tensor of 
the eigenspinors gives the Jacobian matrix of the trans- 
formation from to A; namely it gives the variation of 
the eigenvalues for a small change in the geometry. This 
fact suggests that we can study the map A locally by 
studying the space of the eigenspinor's energy-momenta. 
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By combining ([lC|jl^) and (y_8|) we obtain our main 
result: 

{A", A"} = fd^xfdhj T^-'iix) P^'l{x,y) T^^'^M (20) 



which gives the Poisson bracket of two eigenvalues of the 
Dirac operator in terms of the energy-momentum tensor 
of the two corresponding eigenspinors and of the prop- 
agator of the linearized Einstein equations. The r.h.s. 
does not depend on the gauge chosen for P. 

Finally, if the transformation between the "coordi- 
nates" EI^{x) and the "coordinates" A" is locally invert- 
ible on the phase space T, we can write the symplectic 
form directly in terms of the A"'s as 



n = n„,„ dA" A dA" 



(21) 



where a sum over indices is understood, and where f2„m 
is defined by 



(22) 



Indeed, let dEj^{x) be a (basis) one-form on F, namely the 
infinitesimal difference between two solutions of Einstein 
equations, namely a solution of the Einstein equations 
linearized over E. We have then 

n = Jd^xJ d^y ^lTj{x,y) dElix) A dEi{y) 

= Jd^^J d^y^mn T^n^) dE;,ix)AT"^''j{y) dE^y) 

= ^mn dA" A dA™. (23) 

An explicit evaluation of the matrix would be of 

great interest. 

IV. ACTION 



As shown in 
pressed as 



the gravitational action can be ex- 



s = TT[xm 



(24) 



in natural units fi = G = c = 1. Here x is a smooth 
monotonic function on _R+ such that 



xix) = 1, for .T < 1 — (5, 
x{x) = 0, for X > 1 + 5. 



(25) 



where 5 <<< 1. Equivalently, S is the number of A"'s 
smaller than the Planck mass, which is 1 in natural units. 

The action ( p^ ) approximates the Einstein- Hilbert ac- 
tion with a large cosmological term for "large-scale" met- 
rics with small curvature (with respect to the Planck 
scale). This can be seen as follows. Let N[E] be the 
integer such that A^[i?] is the largest eigenvalue of D[E] 
smaller than the Planck mass Mp = -J— = 1. A moment 
of reflection shows then that we have 



S[E] = N[E] 



(26) 



by definition. For large n, the growth of the eigenvalues 
of the Dirac operator is given by the Weyl formula; 



A" = V- 



(27) 



where V is the volume and we are neglecting factors of 
the order of unity. The next term in this asymptotic 
approximation can be obtained from the fact that the 
Dixmier trace (the logarithmic divergence of the trace) 
of is the Einstein-Hilbert action therefore 



(28) 



(A„)-" = V-^ n-5 + / n 



Under our assumptions on the geometry, at the Planck 
scale we are in asymptotic regime: the first term domi- 
nates over the second, and the remaining terms are neg- 
ligible. Writing the last equation for n = N and using 
(Pq), we have 



1 = V^^ 5- 



Solving for S we obtain 



(29) 



(30) 



which shows that the action ( p^ ) is dominated by the 
Einstein-Hilbert action with a cosmological term. In the 
last equation, we have explicitly reinserted the Planck 
length Lp. 

The presence of the huge Planck-mass cosmological 
term is a bit devastating because the solutions of the 
equations of motions have Planck-scale Ricci scalar, and 
therefore they are all out of the regime for which the 
approximation taken is valid! 

However, the cosmological term can be canceled easily 
by replacing x(x) with xi^) defined by 



X{x) = x{x) - e**x(e2^) 



(31) 



for a small number e. The additional term cancels exactly 
the cosmological term, leaving only the Einstein Hilbert 
action, with corrections which are small for low curvature 
geometries, which, now, are solutions of the theory. In 
fact 

S^TrixiD)) 

= Tt{x{D)) - e^Tr{x{^D)) 

V9R 



V_ 

i^P 



1 

i^P 



_ V 

i^p 



i^p 



jT I VgR 



(32) 
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If we write S directly in terms of the observables A", 
we have the following expression for the action 



(33) 



One cannot vary the A"'s in this expression to obtain 
(approximate) Einstein equations. These equations are 
obtained minimizing ( |3^ ) on the surface X{£), not on the 
entire R°^. In other words, the A"'s are not independent: 
there are relations among them. These relations code the 
complexity of GR. The equations of motion are obtained 
by varying S with respect to the tetrad field. They can 
be written as 



SS _sr' ^A" _ ^ 



dA" 



T"^(x). 
(34) 



Let f{x) — -^x{x). The equations of motion of the the- 
ory are then 



/(A") r"^(x) = 0. 



(35) 



We close by analyzing the content of these equations. 
f{x) is a function that vanishes everywhere except on two 
narrow peaks. A positive peak (width 5 and height 1/5) 
around the Planck mass 1; and a negative peak (width 
(5e and height /5) around the arbitrary large number 
s — ^ » 1. The equation is therefore solved if above 
the Planck mass {n » N), the energy momentum tensor 
scales as 



p{l) T^^(x)=.-^ Pis) T^(^)^(x), 



(36) 



where p{s) is the density of the eigenvalues at the scale 
s and A^*^*^ — s, because in this case the two terms from 
the two peaks cancel. But from (27) we have that the 



density of eigenvalues grows as N^, and that N{s) = 
This yields 



T"^(x) = A" T^^ix) 



(37) 



for n >> N. (Recall that A^ — 1-) In other words: 
the equations of motion for the geometry are solved 
when above the Planck mass the energy-momentum of 
the eigenspinors scales as the mass. 

To understand the meaning of this scaling requirement, 
notice that T"^ is formed by a term linear in the deriva- 
tives of the spinor field and a term independent from 
these, which is a function of (ip, E, dp,E) quadratic in ip. 



(38) 



If we expand the last term around a point of the mani- 
fold with local coordinates x, covariance and dimensional 
analysis require that 



5"^ = cor El +c,Rl + c2REl + 0i^—y (39) 

for some fixed expansion coefhcients cq,ci and C2. Here 
i?^ is the Ricci tensor. [To be convinced that terms of 
this form do appear, consider the following. 

T";, = r^D^r + ■■■ 

= (A")-i ri'rD^D,yj^ + ... 
= (A")-i ^"7V[^M,^^]V'" + ... 

= (A")^i ^w^R^.r + ■■■ 

= (A")-i ^ V7Xf 7J7if ^" + ■ ■ ■ 



Tr I'l^Ri^lJlK 

Ri + ...] 



(40) 



For sufficiently high n, the eigenspinors are locally plane 
waves in local cartesian coordinates: doubling the mass 
just doubles the frequency. If 



A'" = t A" 



(41) 



Then c?^?/''" = t d^ip"-. It follows that in general the 
energy momentum scales as 



(42) 



[ci Rl + C2 REl]+0( — 



For large A" we can disregard the last term, and (|37 
requires that the second square bracket vanishes. Taking 
the trace we have i? = 0, using which we conclude 



r; = o 



(43) 



which are the vacuum Einstein equations. Thus, the vari- 
ation of the (modified) Chamseddine-Connes action im- 
plies a scaling requirement on the high mass eigenspinors 
energy momentum tensors, and this requirement, in turn, 
yields vacuum Einstein equations at low scale. 
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